Origin of Plateau and Species dependence of Laser-Induced High-Energy 

Photoelectron Spectra 
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We analyzed the energy and momentum distributions of laser-induced high-energy photoelectrons 
of alkali and rare gas atoms. For the plateau electrons with energies above 4U P , (U p is the pon- 
deromotive energy), in the tunneling ionization regime, we showed that they originate from the 
backscattering of laser-induced returning electrons. Using the differential elastic scattering cross 
sections between the target ion with free electrons, we explain experimental observations of whether 
the plateau electron spectra is flat or steeply descending, and their dependence on species and laser 
intensity. This quantitative rescattering theory can be used to obtain energy and momentum distri- 
butions of plateau electrons without the need of solving the time-dependent Schrodinger equation, 
but with similar accuracy. 

PACS numbers: 32.80.Rm, 32.80.Fb, 34.50.Rk 
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When atoms or molecules are placed in an intense laser 
pulse, an electron can be released through either a mul- 
tiphoton or a tunneling mechanism. The distinction is 
based on the Keldysh parameter 7 = ^I p /2U P , where I p 
is the ionization energy and U p — Aq/4 is the pondero- 
motive energy with Aq is the peak value of the vector 
potential. In the multiphoton regime (7 > 1), the elec- 
tron spectra exhibit characteristic above-threshold ion- 
ization (ATI) peaks separated by photon energy, with 
yields decreasing monotonically with increasing electron 
energy [Jj]. At higher intensities, in the tunneling region 
(7 < 1), the spectra are notably different. First, the ion- 
ization yield drops steeply from the threshold, but from 
about 3 or 4JJ P onward, the yield flattens out significantly 
until at about \0U P where it drops precipitously again. 
The flattened spectral region from 4-10£/p is called the 
plateau. Similar plateau and cutoff are also well-known in 
the high-order harmonic generation (HHG). Despite this 
canonical description, the HHG and electron spectra in 
the plateau region are not always flat. For electron spec- 
tra, earlier experiments show pronounced enhancement 
in the plateau region in potassium, but not in sodium 
0. Similarly, using nearly identical lasers, clear plateau 
shows up in Xe target, but not in Kr and Ar [3]. Exper- 
imentally the energy spectra of plateau electrons have 
been observed to depend on laser intensities [3| ■ The ori- 
gin of plateau electrons and their dependence on target 
species and laser intensity is explained in this Letter. 

Considerable understanding of laser-induced ATI spec- 
tra from atoms has been achieved since 1990's @, @ 01- 
While the observed spectra in general can be reproduced 
from solving the time-dependent Schrodinger equation 
(TDSE) within the single active electron approximation, 
interpretation of the plateau electrons is based on the 



rescattering model. In this model electrons that are re- 
leased earlier by tunneling ionization can be driven back 
by the laser field to recollide with the target ion. The 
plateau electrons are understood to originate from the 
elastic back scattering of the returning electrons by the 
target ion. However, existing rescattering model does 
not predict how they depend on target species and laser 
intensity. In this Letter we provide such a quantitative 
rescattering theory (QRS) where the energy and momen- 
tum spectra of the plateau electrons can be directly calcu- 
lated using the elastic scattering cross sections between 
the target ion and free electrons, in combination with 
laser- induced returning electron wave packet. We also 
show that the returning electron wave packets depend 
largely on the laser parameters only. Thus the behav- 
ior of the plateau electron spectra is determined entirely 
by the energy and species dependence of the elastic scat- 
tering cross sections between the returning electrons and 
the target ion. 

To understand electron energy spectra, we first exam- 
ine how electron-ion collisions contribute to the photo- 
electron momentum distributions. In Fig. 1 we show a 
typical two-dimensional (2D) electron momentum spec- 
tra for an atomic target calculated by solving the TDSE 
B Hi- (We use linearly polarized electric field for the 
laser pulse along the z axis with the carrier-envelope 
phase set to zero [8|). The horizontal axis is the direc- 
tion of laser's polarization and the vertical axis is any 
direction perpendicular to it (the electron spectra has 
cylindrical symmetry). Due to the short laser pulse, the 
right-hand side (p z > 0) and the left-hand side (p z < 0) 
spectra are not symmetric. In Fig. 1, four semi-circles 
with photoelectron energies of 2, 4, 6 and 10U P , respec- 
tively, are shown. On the right-hand side, three other 
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FIG. 1: Typical 2D electron momentum distributions. Photo- 
electrons with constant energy are given by concentric circles 
from the origin. The three half circles on the right-hand side 
depict momentum surfaces resulting from the elastic scatter- 
ing of returning electrons by the target ion in the presence 
of a laser field. High-energy electrons are obtained only via 
large angles backscattering. 



semi-circles are given, each with its center shifted along 
the p z axis. Each of this circle can be expressed vectori- 
allybyp„ = — A r p z +p r p. Measured from its own center, 
the circle maps out the momentum space of the elastically 
scattered electron with incident momentum p r . The cen- 
ter is shifted since collision occurs in the presence of the 
laser field. If the electron-ion collision occurs at time t r 
when the vector potential is A(t r ) = A r p z , the electron 
will reach the detector gaining an additional momentum 
—A r . (Atomic units are used in the above expressions 
such that the momentum gain in the laser field is di- 
rectly given by the vector potential A r .) In the figure, 
we have chosen the "incident" electrons to be from the 
right, i.e., with negative p z . The direction of p and the 
clastic scattering angles 6 are measured from this "inci- 
dent" direction. This choice would result in high-energy 
or plateau electrons to emerge with positive p z after the 
electron undergoes large-angle backscatterings (9 > 90°). 
Note that electrons scattered into the forward directions 
will emerge with low energies. Clearly, similar circles can 
be drawn for electrons "incident" from the left. 

According to the classical model, electrons that return 
to the ion core at the time when the vector potential is 
at the peak, Aq, it will have maximum kinetic energy 
of 3A7U P = pl/2, i.e., with p r = 1.26 A,. If the elec- 
tron is backscattered by 180°, the emerging photoclcc- 
tron would have a momentum of 2.26^4oj or energy of 
10Z7p. In a recent paper [1 01 ] . we examined the photoelec- 
tron distributions after the p r — 1.26 A electrons have 
been elastically scattered into the backward directions. 
The electron yield along this ring (called BRR, or back 
rescattering ridge) has been shown to be given by 



I{Pv,6v) = S(p r )a(p r ,9). 



(1) 



(DCS), <r(p r ,8), multiplied by a returning electron wave 
packet, represented by S(p r ), with the target ion. The 
validity of Eq. (JTJ) was established in [Icj based on accu- 
rate TDSE results. Its validity has a lso been confirmed 
experimentally in two recent reports 
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This equation is based on the concept of rescattering the- 
ory. It states that photoelectron yields along the BRR 
are proportional to the elastic differential cross section 



To examine the plateau electrons, elastic scattering by 
electrons with returning energy less than 3A7U P must 
be included. These are electrons that return to the par- 
ent ion when A r = A(t r ) is not at the peak. We set 
the returning electron's momentum p r = 1.26 A and 
using the rescattering model, Eq. (fTJ). The validity of 
this model can be checked by comparing the extracted 
a(p r ,9) from the calculated photoelectron momentum 
spectra I(jp v ,6 v ), with those directly calculated from 
electron-ion collisions, as carried out in Morishita et al. 
[l(J. Alternatively, we can first use Eq. (fTJ) at an arbi- 
trary angle, say, for photoelectron with momentum vec- 
tor 10° from the polarization axis, to derive the wave 
packet S(p r ). We then obtain the rescattering "gener- 
ated" photoelectron spectra by using Eq. (fT]). This way 
we have a quantitative rescattering theory which can be 
used to obtain the momentum and energy spectra for 
plateau electrons. We remark later that S(p r ) can also 
be obtained by other approximations. 

In Figs. 2(a,b) electron momentum spectra (p z > 0) 
for sodium and potassium atoms exposed to a 5-cycle 
laser pulse with wavelength of 3200 nm and peak inten- 
sity of 1 TW/cm 2 are presented. On the left-hand half 
of the 2D plot, the results calculated from TDSE are 
shown, and on the right-hand half the same distributions 
obtained from the present QRS model are presented. If 
the two calculations agree well, then there should exist 
a good reflection symmetry in the figure. This is clearly 
the case for each target, confirming the validity of the 
QRS model at the level of electron momentum distribu- 
tions for the plateau electrons. In Figs. 2(c,d) the DCS 
for e-Na + and e-K + collisions are shown for the range 
of momenta p r that contribute to the 4-10U p photoelec- 
trons. As the collision energy increases, the DCS de- 
creases. However, for the same p r , the DCS for e-K + 
is about 100 times larger. Furthermore, the angular de- 
pendence of the DCS in Fig. 2(c,d) are clearly reflected 
in the momentum distributions in Figs. 2(a,b). In other 
words, the elastic scattering cross sections between free 
electrons and target ion can be "read" out directly from 
the ATI electron spectra in the momentum space. This is 
possible because of the other important feature of QRS: 
the returning electron wave packet S(p r ) depends mostly 
on the lasers only. This fact has been noted in the high 
energy region before [13j . 

From the momentum spectra the ATI energy spectra 
can be calculated. Since the right-hand side and left- 
hand side 2D momentum spectra are not symmetric for 
short pulses, we calculate the energy spectra from each 
side separately and the sum of them reproduces the whole 
energy spectra in which the left-hand side makes more 
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FIG. 2: Single ionization of Na and K atoms by a 5-cycle, 3.2 
fj,m infrared laser with peak intensity of 1.0 x 10 12 W/cm 2 
(Up = 0.96 eV). (a,b) Comparison of 2D momentum distri- 
butions above 4[/ p and p z > (i.e., half of the momentum 
space only) obtained from TDSE with those from QRS model. 
Agreement between the two calculations is indicated by the 
good reflection symmetry in each figure; (c,d) Elastic differ- 
ential electron-ion scattering cross sections at large angles; (e) 
Electron energy spectra calculated from TDSE and those from 
the QRS model (dotted lines) above 4U P for each target; (f) 



Experimental electron spectra taken from Gaarde et al. 
In (e,f) the spectra from the two targets are normalized to 
each other at low energies. 



contribution to the lower energy part. In Fig. 2(e), the 
results from TDSE and from QRS are compared. We 
limit the QRS model to electron energy above AU P only. 
Clearly the QRS model reproduces the TDSE results 
well. These calculations can also be compared to the ex- 
perimental data reported in Gaarde et al. |2| , as shown in 
Fig. 2(f). There is a close similarity in the relative yields 
between our calculations and experiment, even though 
our calculations used a 5-cycle (20 fs FWHM) pulse while 
the experiment used somewhat different intensities for a 
pulse with duration of 1.9 ps. The QRS results clearly es- 
tablish that the experimental plateau electron spectra for 
Na and K are due to their elastic scattering cross sections 
in the momentum range of 0.3 to 0.5. Although it has 
long been understood (or speculated) that the behavior 
of plateau electrons are related to elastic scattering cross 
sections [H, 0] , no direct quantitative connection between 
the two has ever been established until now. 

As a side note, we comment that the lower end of the 
plateau electrons has been set at 4U P in this Letter. Ac- 
cording to the classical theory, direct ionization by tun- 
neling (without rescattering) will give maximum electron 




4 6 8 10 12 
Energy/U p 



80 100 120 140 160 180 
Scattering angle (deg) 



FIG. 3: Single ionization of Xe, Kr and Ar atoms by a 5-cycle 
laser pulse. Except for (c), laser's wavelength is 800 nm and 
the peak intensity is I = 1.0 x 10 14 W/cm 2 (A = 0.94). (a) 
Comparison of electron energy spectra from TDSE, and from 
the QRS model above 4{7 P (chain lines). The QRS energy 
spectra for Kr and Xe are calculated using the wave packet 
extracted from Ar, with proper normalization at one energy 
point; (b) Normalized momentum distributions of the return- 
ing electron wave packets extracted from the "left-side" of the 
momentum spectra; (c) Comparison of Xe plateau electron 
spectra at different laser intensities and wavelengths, with 
laser parameters indicated in the figure. Note that the QRS 
model fails at the lowest intensity; (d,e,f) Differential electron- 
ion elastic scattering cross sections used in the QRS model. 



energy of 2U p . Between 2 and 4U P , the direct ionization 
part could still interfere with the rescattering part [1 31 ] . 
thus we chose 4U P as the lower end of the QRS model. 

We have also studied the plateau electrons for Ar, Kr 
and Xe atoms. Experimental data for these systems 
have been reported using lasers of different wavelengths 
and pulse durations since the 1990's [!, H[. In fact, the 
plateau in the high-energy ATI spectra was first observed 
in Xe. In Fig. 3(a) we show the calculated ATI spectra 
vs electron energy in units of U p at the same laser inten- 
sity of 1.0 x 10 14 W/cm 2 for a five-cycle, 800 nm pulse. 
In the case of Xe, the plateau is clearly seen - it re- 
mains almost constant from 5U P to 1QU P , while for Kr 
and Ar, each drops about two to three orders of mag- 
nitude in the same energy range. Since the momentum 
distributions of the returning electron wave packets are 
essentially identical (up to a normalization due to the 
different tunneling ionization rates) for the same laser 
pulse, as shown in Fig. 3(b), the differences in the elec- 
tron energy spectra are attributed to differences in the 
DCS at large angles among the three targets. In Figs. 
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3(d-f) their DCS are shown for the relevant range of p r . 
Note that they cover the same range of magnitude, but 
at large angles, say 160°-180°, the DCS for Xe behaves 
"anomalously" for it increases with increasing energies 
instead of otherwise, as in Ar and Kr. Since higher re- 
turning energies contribute more to the plateau electron 
yields (see Fig. 1), this "anomalous" energy dependence 
explains the flat plateau in Xe. Likewise, the "normal" 
energy dependence of DCS in Ar and Kr can be used to 
explain the steep drop of their electron spectra at high 
energies, see Fig. 3(a). We comment that the QRS re- 
sults shown in Fig. 3(a) are again in good agreement with 
those obtained from TDSE. 

According to the QRS model, for different lasers with 
identical U p , their returning electron wave packets will 
have the same range of kinetic energy, thus similar 
plateau electron momentum and energy distributions. In 
Fig. 3(c), we show the electron spectra for Xe using 800 
nm and 400 nm lasers, respectively, but with the latter 
having four times the peak intensity. The spectra in the 
plateau region indeed agree quite well quantitatively af- 
ter they are normalized to each other at around 10U p . 
Both results are obtained from solving the TDSE. 

We also comment that the QRS model is based on the 
rescattering theory, thus it would fail at lower intensi- 
ties in the multiphoton ionization regime. In Fig. 3(c) 
we show the comparison of results from TDSE and from 
QRS at two lower laser intensities, with Keldysh param- 
eters of 7 = 1.59 and 2.24, respectively. We see clear 
evidence of the failure of the QRS model at 7 = 2.24. 
This sets a rough upper limit for the validity of QRS at 
7 near about 2.0. 

The above examples illustrate that the behavior of 
laser-induced plateau electrons with energies from AU P 
to 10[/p are determined by electron-ion elastic scattering 
cross sections for free electrons with energies between 
\.24U V and 3.17U P (or momentum between 0.79-1.26^4o)- 
Flat plateau is expected when the DCS at large angles 
(close to 180°) increases with increasing kinetic energy 
and when they are highly peaked at large angles, as in 
the case of Xe [Fig. 3(d)]. Such conditions occur often 
in the DCS in low energy electron-ion collisions. Where 
does it occur depends on electron energies and target 
species. Based on the DCS shown in Fig. 3, the plateau 
will become more pronounced for Ar and Kr, but less so 
for Xe, as the laser intensity is decreased. This is consis- 
tent with earlier experimental results 

In summary we studied laser-induced high-energy 
plateau electrons. Together with our previous results 
where the species dependence of HHG was traced to their 



photo-recombination cross sections [13, EH EH , we now 
have established a quantitative rescattering (QRS) the- 
ory for laser-induced ATI electron and HHG spectra in 
the plateau region. Based on the QRS model which is 
valid in the plateau region, there is no need to solve 
TDSE directly. For HHG one only needs to calculate 
photo-recombination cross sections and for plateau elec- 
trons one needs only to calculate elastic electron-ion scat- 
tering cross sections. The returning electron wave pack- 
ets, as shown previously HI, lH, can be extracted from 
a companion target or from the second-order strong-field 
approximation. Conversely, the QRS theory allows one 
to extract electron and photon scattering information 
from laser-induced electron spectra or HHG, respectively. 
As proposed elsewhere [Hj], these cross sections can be 
further used to deduce the structure of the target, thus 
opening up the opportunity of using infrared laser pulses 
for determining the structural change of a dynamic sys- 
tem with temporal resolution of sub-femtoseconds to a 
few femtoseconds. 
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